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1 Introduction 

CLEO, BABAR, and Belle have searched for lepton- number- violation (LNV) in B 
meson decays at e + e _ colliders. In the Standard Model (SM), lepton number L is 
conserved in low-energy collisions and decays, and the lepton flavor numbers for the 
three lepton families are conserved if neutrinos are massless. A number of mechanisms 
for LNV have been proposed including multi-Higgs-boson extensions, leptoquarks and 
Majorana neutrinos. 

A common theme among the experiments is the use of two discriminating variables 
calculated in the center-of-mass frame (CM) (AE = E* B — y/s/2 and mss = m C a,nd = 
Mbc = \Js/A — p* B , where E* B and p* B are the CM energy and momentum of the 
reconstructed B meson candidate), and event-shape discriminants. The signal events 
peak at zero for AE, at approximately the B meson mass for rngg, and the event- 
shape is spherical; for the backgrounds, AE and mgs are smoothly varying, while 
the event-shape is more jet-like. The major backgrounds are continuum production 
of quark pairs e + e~ — > qq (q = u, d, s and c), semileptonic B decays which produce a 
lepton, and other BB decays where a hadron is mis-identified as a lepton. 



2 Recent results from CLEO, BABAR, and Belle 

CLEO reused the techniques for their search for B — > X s £ + £~ to look for B + — > 
h~e + e + , h~e + fi + and h~fi + fi + where h~ = K~ , K*~ ,tt~ and p~ [TJ, based on 9.6 
million BB events. There are three main sources of background: B decays of the 
type B — > XJ/ip and B — > Xip(2S); other B decays, with two apparent leptons 
(either real leptons or hadrons misidentified as leptons); and continuum processes 
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with two apparent leptons. The backgrounds from J/ip and ip{2S) were severe in the 
searches for B — > X£ + £~ as the decays provided two opposite sign leptons but here 
they contribute less than 0.1 event per decay mode after careful lepton identification. 

An unbinned maximum likelihood (ML) method is used to discriminate between 
signal events and the remaining two background sources. The variables are: the 
missing event energy, -E m i SS = 2E heiira — -^det, where J2E^ e t denotes the sum of 
energies of all the detected particles in the event; a Fisher discriminant J 7 based on the 
event-shape; M cand ; and AE. Loose cuts of l-Ewsl < 2 GeV, 5.2 < M cand < 5.3 GeV/c 2 
and \ AE\ < 0.25 GeV are applied. The branching fraction for the signal and the yields 
for the two backgrounds are free parameters in the ML. No evidence for the decays is 
found and all modes have a statistical significance of less than 1.2 standard deviations. 
90% confidence level (CL) upper limits (UL) are placed on the branching fractions as 
shown in Table [TJ The upper limits range from 1.0 to 8.3 x 10 -6 . 
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Table 1: CLEO results with the statistical significance S and 90% confidence level 
upper limits UL (xl0~ 6 ) on the branching fraction (including systematic error) [TJ. 

BABAR has repeated the CLEO search for the four modes B + — > h~ p + p + and 
B + — > h~e + e + , where h~ = K~ or tt~ [2] and they also reuse the techniques developed 
for their B — > K^£ + £~ analyses, using a data sample of 471 ± 3 million BB pairs. 
The key difference with CLEO is the ML function used. The two main backgrounds 
from continuum events and BB decays are suppressed through the use of a ML ratio 
1Z constructed from boosted decision tree discriminants (BDTs). Before the ML fit, 
a selection on TZ retains 85% of the simulated signal events while rejecting more than 
95% of the background. The selection efficiency for simulated signal is 13%-48%. 

The signal branching fraction and background yields are extracted from the data 
with an unbinned ML using mss and TZ. The signal rriEs distributions are taken from 
data using a Gaussian shape unique to each final state, with the mean and width 
determined from fits to the analogous final states in the B + — > J/ip(-^r l + l~)h + 
events from the data. No significant yields are observed and the results of the ML 
fits to the data are summarized in Table [2j 

Belle have performed the first searches for the decays B + — > D~e + e + , D~e + p + 
and D~p + p + using a data sample containing 772 million BB pairs [3]. As the CKM 
matrix element \V c b\ > \V u b\, these decays could be enhanced by an order of magnitude 
compared to B + — > h~l + l + . Belle look for an energetic same-sign dilepton and 
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Table 2: BABAR results, showing the total events in the sample, signal yield and its 
statistical uncertainty, significance S, branching fraction B, and 90% CL branching 
fraction upper limit Bjjl 0- 

combine it with a D candidate requiring a proper charge combination for the dilepton. 
The same-sign lepton pair must have a total energy in the T(4S) center-of-mass (CM) 
frame greater than 1.3 GeV. More than 95% of events have only one same-sign lepton 
pair. 

Candidate D~ mesons are reconstructed in the D~ — > K^ix-ix^ decay. The three 
tracks from the D~ candidate are fit to a common vertex and are required to have a 
7^ + 7r~7r~ invariant mass {Mk-k-k) within approximately ±10 MeV/c 2 from the nominal 
D~ mass. The Mk-k-k distribution is fit to two Gaussian functions with a common 
mean in a mass window ±3 times the width of the narrower Gaussian component. 
The average multiplicity of D~ candidates is 1.3 per event. 

The continuum background is rejected using a Fisher discriminant J 7 based on 
event-shape variables, and cos6b, the cosine of the polar angle of the B candidate 
flight direction in the CM frame. Semileptonic decays such as B — > D~£ + i>fX have 
missing energy due to the undetected neutrino and the two reconstructed leptons do 
not come from the same vertex. These events can be rejected using the missing energy 
E m i SS and Sz, the separation between the impact parameter of the two leptons in the 
beam direction. The four variables, J 7 , cos 9b, E miss and Sz, are combined together 
into a single likelihood ratio, R s . The optimal requirement on R s is determined by 
maximizing the figure of merit, e s /(a/2 + y/N^), where e s is the MC signal efficiency, 
iVb is the number of expected background events in the signal region, and a is set to 
zero. 

After applying the R s requirements, 5, 23 and 40 events remain in the background 
region for the e + e + , e + /z + and modes, respectively. The signal efficiencies are 

evaluated to be 1.2% - 1.9%. The expected numbers of background events in the 
signal region are 0.18, 0.83 and 1.44 events for the e + e + , e + fi + and modes, 
respectively. Figure Q] shows the M^ C -AE distributions of the selected events. No 
events are observed in the signal region. Table |3] summarises the upper limits achieved. 

BABAR has performed the first measurement of the branching fraction upper lim- 
its for the baryon- and lepton- number violating decays B° — > A~^£~, B~ — > Ai~ and 
B~ — > At~ [4j with 471 ± 3 million BB pairs. The and A candidates are recon- 
structed through the decay modes A+ — > r pK~it + and A — > p7v~ , respectively. The 
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Figure 1: The M^ C -AE distributions of D e + e + (left), D e + /i + (middle) and 
D~fi + fi + (right) final states in data. The (red) boxes indicate the signal regions. 
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Table 3: Belle results for the B + — > D~£ + £ + search; e is the signal reconstruction effi- 
ciency; iV bs is the number of events in the signal region; is the expected number 
of background events, and UL is the branching fraction 90% CL upper limit [3]. 

final state tracks for both the and A decays are constrained to a common spatial 
vertex and their invariant mass is constrained to the A^ or A mass. The candidates 
are required to be within ±15 MeV/c 2 of the nominal A+ mass and ±4 MeV/c 2 of 
the nominal A mass. B meson candidates are formed by combining baryon candidate 
with a or e~ and constraining them to a common point. Bremsstrahlung energy 
recovery is performed for electrons. Background from e + e~ — > e + e~7 events in the 
Al channel are eleminated by requiring more than four charged tracks in the event. 

Candidate selection is optimized using a figure of merit (see above) with a = 5. A 
neural net (NN) is used to provide further discrimination between signal and back- 
ground such that about 90% of the signal is retained and about 50% of the background 
is rejected. The remaining background after selection for the A+£~ modes is com- 
posed of roughly equal amounts of BB and continuum events, while the background 
for the Ai modes is almost entirely continuum. 

The variables AE and 7Bes are used in the ML fit to the two Al modes; the A+£~ 
decay has more background and the NN is used as a third discriminating variable. 
No significant signal is observed and an upper limit is calculated for the branching 
fraction for each decay mode, as shown in Table HI 
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Table 4: Summary of number of candidates (iV can d), branching fraction central value 
(B), signal efficiency (e), and branching fraction 90% CL UL (B 90 %) [I]. 
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Figure 2: Branching fraction upper limits for B meson LNV decays from CLEO pQ, 
BABAR [21 S], Belle [3], and LHCb (5]. 
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